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ABSTRACT 

We measure the correlation between the arrival directions of the highest energy cosmic rays 
detected by the Pierre Auger Observatory with the position of the galaxies in the H I Parkes All 
Sky Survey (HIPASS) catalogue, weighted for their H I flux and Auger exposure. The use of 
this absorption-free catalogue, complete also along the galactic plane, allows us to use all the 
Auger events. The correlation is significant, being 86.2% for the entire sample of H I galaxies, 
and becoming 99% when considering the richest galaxies in HI content, or 98% with those 
lying between 40-55 Mpc. We interpret this result as the evidence that spiral galaxies are the 
hosts of the producers of UHECR and we briefly discuss classical (i.e energetic and distant) 
long Gamma Ray Burst (GRBs), short GRBs, as well as newly born or late flaring magnetars 
as possible sources of the Auger events. With the caveat that these events are still very few, and 
that the theoretical uncertainties are conspicuous, we found that newly born magnetars are the 
best candidates. If so, they could also be associated with sub-energetic, spectrally soft, nearby, 
long GRBs. We finally discuss why there is a clustering of Auger events in the direction on 
the radio-galaxy Cen A and an absence of events in the direction of the radio-galaxy M87. 

Key words: cosmic rays - galaxies: gamma-rays: bursts - galaxies: statistics — radio lines: 
galaxies 



1 INTRODUCTION 

The origin of ultra-high energy cosmic rays (UHECR), exceeding 
10 EeV (1 EeV=10 18 eV) has been a mystery for decades, but the 
recent findings of the large area detectors, such as AGASA (Ohoka 
et al. 1997), HIRes (Abu-Zayyad et al. 2000), and especially the 
Pierre Auger Southern Observatory (Abraham et al. 2004), began 
to disclose crucial clues about the association of the highest energy 
events with cosmic sources. The Auger collaboration (Abraham et 
al. 2007) found a positive correlation between the arrival directions 
of UHECR with energies grater than 57 EeV and nearby AGNs (in 
the optical catalogue of Veron-Cetty & Veron 2006). Although this 
result has not been confirmed by HIRes (Abbasi et al. 2008) and it 
has been criticised by Gorbunov et al. (2008), it received an impor- 
tant confirmation by George et al. (2008), who considered a com- 
plete sample of nearby hard X-ray emitting AGNs detected by the 
BAT instrument onboard Swift. This sample is much less affected 
by absorption than any optical sample although, to identify as such 
an AGN, one relies on optical identification. Moreover, George et 
al. (2008) found a correlation not simply with the AGN locations, 
but by weighting them for the X-ray flux and the Auger exposure. 

This association, if real, is surprising, since the large ma- 
jority of the correlating AGNs are radio-quiet, a class of objects 
not showing, in their electromagnetic spectrum, any sign of non- 
thermal high energy emission: no radio-quiet AGN was detected by 
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the EGRET instrument onboard the Compton Gamma Ray Obser- 
vatory (Hartman et al. 1999). Therefore they must accelerate par- 
ticles (protons, nuclei, and presumably their accompanying elec- 
trons) to ultra-high energies without any noticeable radiative emis- 
sion from these very same particles. Radio-loud AGNs, instead, 
together with Gamma Ray Bursts (of both the long and short cat- 
egory) do show high energy non-thermal emission, and have been 
considered for a long time better candidates as UHECR sources 
(Vietri 1995; Waxman 1995; Milgrom & Usov 1995; Wang, Raz- 
zaque & Meszaros 2008; Murase et al. 2008; Torres & Anchor- 
doqui 2004 and Dermer 2007 for reviews, and Nagar & Matulich 
2008 and Moskalenko et al. 2008 for the possible association of the 
AUGER events with radio-loud AGNs). Note also that some short 
GRBs could be due to the giant flares of highly magnetised neu- 
tron stars ("magnetars", as the 27 Dec 2004 event from 1806-20; 
Borkowski et al. 2004; Hurley et al. 2005; Terasawa et al. 2005), 
and that, at birth, a fastly spinning magnetar can be much more 
energetic than when, later, it produces giant flares (Arons 2003). 

The possibility that GRBs and magnetars are the sites of pro- 
duction of UHECR would directly imply the direct association of 
these events with (normal) galaxies. In this case the found associ- 
ation of UHECR with nearby AGNs might then be due to the fact 
that local AGNs just trace the distribution of galaxies. The aim of 
the present paper is to test this possibility directly correlating the 
locations of the ultra-high energy Auger events with a well de- 
fined, complete, and possibly absorption-free sample of galaxies. 
For this purpose we use the sample of HI emitting galaxies, com- 
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piled using the Parkes 64- m radio telescope (Barnes et al. 2001; 
Staveley-Smith et al. 1996), which is conveniently located in the 
south hemisphere, as the Auger observatory. The entire sample cov- 
ers the portion of the sky visible by Auger, making it possible to 
use, for the correlation analysis, all the 27 UHECR events with en- 
ergies larger than 57 EeV detected by Auger, without excluding the 
galactic plane, as is instead necessary when dealing with AGNs or 
optically selected galaxies. Note that the presence of neutral hydro- 
gen strongly favours spirals (or, more generally, gas-rich galaxies) 
with respect to elliptical galaxies. 

We use a cosmology with ho = Qa = 0.7 and S7m = 0.3. 



2 DATA 

2.1 UHECR events 

The Auger Observatory (Abraham et al. 2004, 2008), operating in 
Argentina since 2004, is located at latitude —35.2° and it has a 
maximum zenith angle acceptance of 60°. The relative exposure is 
independent of the energy of the detected events and it is a nearly 
uniform in right ascension. The dependence on declination is given 
by Sommers (2001). The Observatory can detect Cosmic Rays from 
sources with declination 8 < 24.8°. 

The available Auger list of UHECR events (Abraham et al. 
2008) comprises 27 events with energies in excess of 5.7 x 10 19 eV 
from an integrated exposure of 9000 km 2 sr yr. The event arrival 
directions are determined with an angular resolution of better than 
1°. However, magnetic fields of unknown strength will deflect 
charged particles on their trajectories through space. The advan- 
tage of studying the highest energy events is that this deflection is 
minimised, but it can still be up to ~ 10° in the Galactic field. 
The 27 UHECR detected by Auger are distributed in the range 
8 e [-61,9.6] (or at galactic latitudes b e [-78.6,54.1] - open 
circles in Fig. Q] and Fig. [2}. 

2.2 HIPASS catalogue 

We compare the arrival directions of Auger UHECRs with the loca- 
tions of sources of the H I Parkes All-Sky Survey (HIPASS - Meyer 
et al. 2004). This is a blind survey of sources in H I covering the full 
southern sky at 8 < 25° which is the same sky area accessed by 
the Auger Observatory. The full catalogue is composed by a list of 
4315 sources at 8 < 2° (HICAT - Meyer et al. 2004; Zwaan et al., 
2004) and by its extension to the northern sky up to 8 — 25° (NHI- 
CAT - Wong et al. 2006) which includes 1002 sources. All sources 
are shown in Fig. Q] with the 27 UHECR detected by Auger. 

The HICAT and NHICAT have different level of complete- 
ness. To have a catalogue complete in flux at the 95% level, we cut 
the HICAT at S int > 7.4 Jy km s _1 and the NHICAT at S int > 15 
Jy km s" 1 as discussed in Zwaan et al. (2004) and Wong et al. 
(2006). Sixxt represents the total HI line flux. For the purposes of 
this paper we also considered the H I sources within 100 Mpc which 
is the maximum distance at which UHECRs of E > 57 EeV can 
survive the GZK suppression effect (see e.g. Harari et al. 2006). We 
call this sample 95HIPASS: it contains 2414 sources from the HI- 
CAT and 290 sources from the NHICAT for a total of 2704 sources 
and covers the entire sky at 8 < 25°. We will also consider the 
southern sky sample alone which is more complete and can be cut 
at 99% completeness for Sint > 9-4 Jy km s -1 (also by consider- 
ing sources at < 100 Mpc). This sample contains 1946 sources and 
is called 99HICAT. 




Figure 1. H I Parkes All-Sky Survey (HIPASS) galaxies (solid filled sym- 
bols) and the 27 AUGER UHECRs (black open circles) in galactic coordi- 
nates. Blue levels corresponds to the AUGER relative exposure. The south 
HI catalogue (HICAT - 4315 sources) and the north extension (NHICAT 
- 1002 sources) are shown by red filled circles and orange stars, respec- 
tively. The green lines are 20° X 20° boxes centred on the positions of the 
radio-galaxies M87 (near the north Galactic pole) and Cen A. 
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Figure 2. HIPASS galaxy H I flux, in galactic coordinates. Blue levels cor- 
responds to the integrated flux (in bins of 2° X 2° and units of ly km s ) 
of H I emission, multiplied by the relative AUGER exposure. Red circles 
are the locations of the 27 AUGER UHECRs above 57 EeV. 

3 ANALYSIS 

To quantify the possible correlation between UHECR Auger events 
and the distribution of H I local galaxies we use the method adopted 
by George et al. (2008). In order to quantify the probability that 
two sets of sources are drawn from the same parent population 
of objects we perform the two-dimensional generalisation of the 
Kolmogorov-Smirnov (K-S) test (Peacock 1983) proposed by 
Fasano & Franceschini (1987). 

In our case the test is used to compare two data samples, i.e. 
the UHECR and the H I galaxies. This test can then measure either 
if UHECRs have a galaxy counterpart, and, viceversa, if a concen- 
tration of galaxies has an UHECR counterpart. The test relies on the 
statistic D, also used for the unidimensional K-S test, which repre- 
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sents the maximum difference between the cumulative distributions 
of the two data samples. For each UHECR data point j we compute 
a set of four numbers dj,i (z=[l,4]) defined as the difference of the 
relative fraction of UHECR and H I galaxies found in the four nat- 
ural quadrants defined around point j. Hence, D — max(dj t i) for 
all the data points considered. Defining Z n — D^/n, the strength 
of the correlation between two catalogues is the integral probability 
distribution P(D v / n > observed), where n = N1N2/IN1+N2), 
and iVi and N2 are the number of data points in the two sets. This 
measurement can be used to determine the similarity of sets of po- 
sitions on the sky. 

The probability can be computed analytically for large data 
sets (n >80 - Fasano & Franceschini 1987). In our case, having 
only 27 UHECR, we have to rely on Monte Carlo simulations. We 
generate a large set of random UHECR events according to the rel- 
ative Auger exposure. For each synthetic UHECR sample we com- 
pute Z n by correlating it with the catalogue of H I galaxies. The 
probability of the observed Z n is given by the number of times we 
find a value of Z n larger than the observed one. This is the probabil- 
ity that the correlation between the (real) UHECR sample and the 
H I galaxies is not by chance. Large (low) values of the probability 
indicate a good (poor) correlation between the Auger UHECRs and 
the given H I galaxy sample. 

As noted by George et al. (2008), the two-dimensional K- 
S test can be performed with the number of data points or with 
the flux of the sources in the comparison sample. In our case D 
represents the maximum difference between the number of UHE- 
CRs and that of the sum of the galaxies weighted for their flux 
and for the the relative Auger exposure. The advantage of using the 
weighted flux of the sources is that it accounts for their distance. 
George et al. (2008) found that the UHECRs are more correlated 
with the weighted flux of Swift AGN than with with their position. 
In Fig. [2] we show the map of the flux of the HIPASS catalogue 
weighted for the Auger relative exposure. 



4 RESULTS 

We found that with the 95HIPASS catalogue (2704 HI sources 
complete in flux at 95%) the probability that UHECRs are cor- 
related with HI galaxies is 71.6% by using the weighted flux of 
the HI sources. Considering the more complete 99HICAT (1946 
HI sources complete in flux at 99%) distributed within 100 Mpc 
and the 25 UHECRs distributed in the same sky region we find 
a larger flux-weighted probability of 87.8%. This probability is 
slightly smaller than found with local AGN by George et al. (2008). 

However, having a large sample of H I galaxies we can study 
if the correlation probability changes by considering different sub- 
samples of galaxies selected according to their distance or lumi- 
nosity. We have considered 4 bins of distance with an equal num- 
ber of sources (~500) per bin. The correlation probability shows a 
maximum of 95% (97.8% for the 99HICAT) for sources distributed 
between 37.8 and 55 Mpc. We show these results in Fig.[3](open 
circles and stars in the bottom panel). 

Similarly we defined four equally populated luminosity bins, 
or, equivalently, four bins of HI mass content, since we can use 
M/Mq = 2.36 x 10 5 D^ pc Sint to estimate the HI mass (here 
Sint is measured in [Jy km/s]). We find that the probability (left 
panel in Fig. [3} is maximised by the most HI luminous or massive 
(in HI) sources (98% and 99% for the 95HIPASS and 99HICAT 
sample, respectively, for M > 1.1 x 10 10 M Q ). 

Selecting those HI galaxies located within two 20° x 20° 



boxes centred on the radio-galaxies Cen A and M87 (green boxes 
in Fig. [TJ, we can show where they lie in the luminosity-distance 
plane in Fig. [3] (orange and green dots, respectively). While there 
is no clustering of points at the distances of Cen A and Virgo, we 
can see that H I galaxies in the direction of Cen A do cluster at dis- 
tances of 40-50 Mpc, where the Centaurus cluster is. This could 
explain why some UHECR events appear to be associated with the 
radio-galaxy Cen A, and none with M87: beyond Cen A there is 
the Centaurus cluster, richer of H I emitting spirals than the Virgo 
cluster. The ratio of the integrated HI fluxes from the two 20° x 20° 
boxes (Virgo/Cen A) is 5.9. To this, we have to multiply by another 
factor 3 for the lower Auger exposure in the direction of Virgo. 

The sample has too few galaxies beyond 100 Mpc to test the 
GZK effect (that would be revealed by finding no correlation for 
these galaxies). 



5 DISCUSSION 

The 27 Auger events above 57 EeV, with a total exposure of 
9000 km 2 sr yr correspond to an integrated flux, in CGS units: 



F A (E > 57 EeV) 



1 1 , 1A-11 -2 -1 

1.1 x 10 erg cm s 



(1) 



This flux is smaller than the electromagnetic flux that we receive 
from nearby radio-quiet AGNs in hard X-rays (see e.g. Tueller et 
al. 2008). We now compare this flux with the expected flux of other 
candidate sources. We will consider flaring or bursting sources, that 
is impulsive events, but the spreading of the arrival times of UHE- 
CRs from a source located at a distance D, At ~ DO 2 /2c, due to 
even tiny magnetic deflections, ensure that we can treat all candi- 
date sources as continuous. We will estimate the predicted flux in 
two different ways. 

First, assume that a class of sources is characterised by a pulse 
of emission of UHECRs, of average total energy < E > . Assume 
also that these events occur at a rate R per galaxy, per year, and con- 
sider those events occurring within the GZK radius D c . We have: 



F =<E> 



R N e (D < D c ) 



3.15 x 10 7 4vr(a_D c ) 2 



(2) 



where 3.15 x 10 7 is the number of seconds in one year and 
N g (D < D c ) is the number of galaxies within D c of L, lumi- 
nosity. The average distance of the sources is aD c (a = 3/4 for 
sources homogeneously distributed). Setting the mean local galaxy 
density n g = 7V g /(47r = 10~ 2 n g ,_ 2 Mpc" 3 we have: 

-Dc.100 



F ~ 1.2 X 10"°' < E > Rn R 



-2 -1 
erg cm s 



(3) 



where D c — 100D c ,ioo Mpc. 

The second estimate on the predicted UHECR flux uses the 
electromagnetic flux as a proxy. Assume that we detect, for a typ- 
ical member of a class of sources, an average fluence < T >, 
and that there are TV events per year. If a fraction r\ of these events 
comes from sources within D c , we have 

%]<T > N 



F = 



3.15 x 10 7 



(4) 



This estimate is more appropriate when dealing with sources, such 
as long and short GRBs, whose fluences and occurrences are 
known, while Eq.[3]is more appropriate when dealing with possible 
sources of unknown electromagnetic output, but predicted energet- 
ics and rates, such as newborn magnetars (Arons 2003) or giant 
flares from old magnetars. 

Let us consider the above classes of sources in 
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Figure 3. Top left panel: the HI mass (in solar masses) of our galaxies as a function of their distance. The HI mass is proportional to the HI luminosity, and 
is found using M = 2.36 x 10 5 D^ pc 5 int , where Sjnt is the integrated flux in [Jy km/s]. Black empty circles are those galaxies forming a complete, flux 
limited, sample. Orange and green filled circles are galaxies in the 20° X 20° boxes centred on the position of the radio-galaxies Cen A and M87, respectively. 
The distances of these two radio-galaxies are marked by an arrow. At the distance of Cen A there is almost no H I emitting galaxy, and no concentration is 
seen at the distance of M87 (and the Virgo cluster). The H I galaxies lying in the same region of the sky as Cen A show a concentration for distances 40-50 
Mpc, where the Centaurus cluster of galaxies is. No concentration in distance is seen for H I galaxies in the direction of Virgo. The bottom left panel shows 
the significance of the 2D K-S test using galaxies in different bin of distances (circles: South sample HICAT; stars: South+North sample HICAT+NHICAT). 
The top right panel shows the significance of the correlations for different bins of H I content (or, equivalently for different luminosity bins). 



turn, starting from short GRBs. In the BATSE catalog 
(cossc.gsfc.nasa.gov/docs/cgro/batse/BATSE_Ctlg/flux.html) 
we have 490 short GRBs of total fluence 5.5 x 10~ 4 erg cm" 2 
in 9 years of operation. Tanvir et al. (2005) correlated these short 
GRBs with local optically selected galaxies finding that a fraction 
between 5 and 25% of BATSE short GRBs might be nearby, i.e. at 
2 < 0.025, corresponding to 109 Mpc. Considering that BATSE 
saw half of the sky and setting r\ — 0.1 we have an average flux 
of 3.9 x 10" 13 erg cm" 2 s" 1 . Then, if the UHECR flux is similar 
to the electromagnetic one, short GRBs do not match the required 
flux. 

Classical long GRBs (namely, energetic GRBs at z>l) in the 
BATSE sample have a total fluence of 0.024 erg cm" 2 (for the 
listed 1490 long GRBs in the BATSE catalog), corresponding to 
an average (all sky) flux of 1.7 x 10 -10 erg cm -2 s _1 , larger than 
the one given by Eq. Q] However, for long BATSE GRBs, r\ must 
be much smaller than 0.1, as directly suggested by the paucity of 
nearby events, and by the lack of correlation with nearby galaxies 
and clusters (Ghirlanda et al. 2006). While we cannot dismiss them 
as sources of UHECRs, it seems likely that classical BATSE bursts 
are too distant (but see below). 

Consider now giant flares from relatively "old" magnetars. 
The giant flare from SGR 1806-20 of Dec 27, 2004 emitted an 
energy E ~ 10 46 erg in less than a second. The radio afterglow 
convincingly demonstrated the formation of a (at least mildly) rel- 



ativistic fireball. With the current hard X-ray instruments, such 
flares can be detected up to ~30-40 Mpc (Hurley et al. 2005). Eq. 
[3] with D cA oo = 0.3, and < E >= 10 46 erg, would require R ~ 1 
event per galaxy per year, while an approximate limit to the rate is 
R < 1/30 yr" 1 (see e.g. Lazzati et al. 2005). 

Finally, consider fastly spinning newly born magnetars, whose 
rotational energy can exceed 10 52 erg, with a rate of R = 10 -4 
events per galaxy per year (Arons 2003). With the estimated galaxy 
density (n g ,_2 ~ 0.7 with L ~ L*; Blanton et al. 2001) there 
should be 1 event per year within 100 Mpc. If each magnetar pro- 
duces 10 50 erg in UHECRs, then this class of sources can be the 
progenitor of the Auger events (Eq.[3]l. This is independent of col- 
limation, since the reduced rate of events pointing at us is com- 
pensated by an increase of the apparent energetics. But if an equal 
amount of energy is released in electromagnetic form, at energies 
detectable by BATSE, then they should be a significant fraction 
of all BATSE GRBs. Since the birth of a magnetar should be ac- 
companied by a supernova, these events should be associated with 
long, rather than short GRBs, for which no associated supernova 
has been seen. If the radiative output is isotropic, they will all be 
nearby, sub-energetic, GRBs. 

The required fluence of these sub-energetic nearby long 
GRBs, to match the UHECRs flux, should be 

< T > ~ 3.15 x 10~ 4 eca FA '~ n erg cm" 2 (5) 
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where ecR is the ratio of the emitted energy in radiation and 
UHECR. If r\ ~ ecR ~ 1, these events constitute a sizeable frac- 
tion of the total fluence of all long BATSE GRBs in one year (which 
is F ~ 0.024/9 ~ 2.7 x 10 -3 erg cm" 2 ). 

Since we know that the large majority of long GRBs are not 
nearby, newly born magnetars should not constitute conspicuous 
events in hard X-rays. Their fluence must be mostly emitted in an- 
other energy range. GRB 060218 (Campana et al. 2006) with an 
energy of a few x 10 49 erg, at a distance of 145 Mpc, could be 
one of these events, and Soderberg et al. (2006) and Toma et al. 
(2007) already suggested that this GRB was powered by a newly 
born magnetar. The spectrum of its prompt emission peaked at ~5 
keV, i.e. its fluence in relatively soft X-rays exceeded the 15-150 
keV fluence. It was also very long, slowly rising, and would not 
have been detected by BATSE. Soderberg et al. (2006) pointed out 
that these sub-energetic long GRBs should not be strongly beamed 
(not to exceed the rate of SN Ib,c), and should occur at a rate of 
230iqo° Gpc -3 yr _1 , corresponding to R « 10~ 5 events per L, 
galaxy per year, about ten times larger than for classical long GRBs 
whose radiation is collimated into 1% of the sky. According to this 
rate, Eq. [3] would then demand < E >~ 6 x 10 50 erg in UHECRs 
to match the observed flux. 



6 CONCLUSION 

We have correlated the cosmic rays with E > 57 EeV detected 
by the Auger Observatory with a complete, absorption-free sample 
of H I selected galaxies. We found a significant correlation when 
correlating the Hi flux of galaxies of our sample. 

When considering the largest 95HIPASS catalogue and the 27 
UHECRs we find a weak correlation (probability of 72%), while 
a larger significance (87.8%) is reached if we consider the most 
complete 99HICAT sample of galaxies (though with 25 UHECRs). 
These probabilities are maximised by cutting the H I sample in dis- 
tance or luminosity bins: it becomes 99% when considering the 500 
most luminous (or most HI massive) galaxies (1/4 of the sample), 
and 98% when considering the 500 galaxies lying between 38 and 
54 Mpc, where the Centaurus cluster of galaxies is. 

Thus there is the possibility that the UHECRs coming from the 
direction of Cen A are instead coming from the more distant Cen- 
taurus cluster. Galaxies of this cluster are richer in H I than Virgo 
galaxies, explaining why there is no UHECR event from the direc- 
tion of Virgo. 

This sample is formed by H I emitting galaxies, therefore it is 
biased against ellipticals. The found correlation with these galaxies, 
per se, is not disproving the found correlation with AGNs (Abra- 
ham et al. 2007, 2008; George et al. 2008), since they also trace 
the local distribution of matter, as spiral galaxies do. On the other 
hand, it opens up the possibility, on equal foot, that UHECRs are 
produced by GRBs or newly born magnetars (see also Singh et al. 
2004 who used AGASA events). With the caveat that it is prema- 
ture, with so few events and big theoretical uncertainties, to draw 
strong conclusions, we have pointed out that although classical (i.e. 
energetic) long GRBs and short GRBs have difficulties in produc- 
ing the required UHECR flux, newly bom magnetars can. If so, 
they could also be a subclass of long GRBs, possibly sub-energetic 
and relatively nearby, powered by fastly spinning, newborn magne- 
tars. The future increased statistics of UHECRs arrival directions 
will help to discriminate among the different proposed progenitors, 
especially if there will be (or not) an excess of events close to the 
radio core and/or lobes of Cen A. 
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